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ABSTRACT 

The nature of the accreting companion to Mira — the prototypical pulsating asymptotic giant branch 
star — has been a matter of debate for more than 25 years. Here we use a quantitative analysis of the 
rapid optical brightness variations from this companion, Mira B, which we observed with the Nickel 
telescope at Lick Observatory, to show that it is a white dwarf (WD). The amplitude of aperiodic 
optical variations on time scales of minutes to tens of minutes (« 0.2 mag) is consistent with that 
of accreting WDs in cataclysmic variables on these same time scales. It is significantly greater than 
that expected from an accreting main-sequence star. With Mira B identified as a WD, its ultraviolet 
(UV) and optical luminosities, along with constraints on the WD effective temperature from the UV, 
indicate that it accretes at ~ 10 -10 Mq yr _1 . We do not find any evidence that the accretion rate is 
higher than predicted by Bondi-Hoyle theory. The accretion rate is high enough, however, to explain 
the weak X-ray emission, since the accretion-disk boundary layer around a low-mass WD accreting 
at this rate is likely to be optically thick and therefore to emit primarily in the far or extreme UV. 
Furthermore, the finding that Mira B is a WD means that it has experienced, and will continue to 
experience nova explosions, roughly every 10 6 years. It also highlights the similarity between Mira AB 
and other jet-producing symbiotic binaries such as R Aquarii, CH Cygni, and MWC 560, and therefore 
raises the possibility that Mira B launched the recently discovered bipolar streams from this system. 
Subject headings: stars: individual (HD 14386, VZ Ceti) — binaries: symbiotic — novae, cataclysmic 
variables — accretion, accretion disks — white dwarfs 
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1. INTRODUCTION 

Variations in the strength of emission from the accre- 
tion region around a compact object can reveal funda- 
mental properties of the compact object and its accre- 
tion structure. For example, the frequencies associated 
with features in the power density spectra (PDS) of these 
brightness variations are linked to physical ti me scales 
at different location s in the accretion disk (e.g.. iMauchd 
120021: IWarnerl l200i IDone. Gierlihski. fc Kubo"tar i2007h 
Such PDS features include quasi-periodic oscillations, 
"breaks" where the slope of the PDS changes, and broad- 
band com ponents that can be fitted with Lor entzian 
functions (Psaltis, Bcll oni. k, van der Klisl [19991 ). Em- 
pirically, the locations of these features depend on the 
type of object; in X-ray binaries, the highest frequency at 
which significant variability exist s is greater for neutron- 
star than black- hole accretors (jSunvaev fc Revnivtsevl 
2000), consistent with the higher Keplerian frequency 
at the inner edge of the accretion disk around a neu- 
tron star. Even an apparently featureless PDS may 
hold the imprint of the underlying accreting object, 
as is evidenced by the correlation between the nor- 
malization of the PDS at a particular frequency and 
the mass of the black hole in active galaxies compared 
to X -ray binaries dNikolaiuk. Papadakis. fe Czernvll2004l : 



iGierlihski. Nikola ink, fe CzernvH2008i r 

An important case where accretion-related variations 
can reveal the nature of an accreting object is the 
wide, interacting binary-star system Mira AB (o Ceti). 
Mira AB is a 'symbiotic-like' or 'weakly symbioti c' sys- 
tem at a distance of 107 pc (jKnapp et al.1 [2003) . Ac- 



cording to the preliminary orbit of iPrieur et al.l ([2002) , 
the orbital period is 497.9 years, and the inclination is 
112.0° (i.e., we see the binary nearly edge-on). The donor 
star is the prototype of the Mira-type class of p ulsating 
AGB stars, with a pulsation period of 332 days (Hofficit 
|1997|) . A more compact companion, whose nature has 
been controversial, accretes from the wind of the AGB 
star. Although mo st indicators suggest that the accre- 
tor is a WD (e.g., IWarnerl Il97l iReimers fe Cassatellal 
|l985j), several authors have argued that the low X- 
ray luminosity implies it is a main-sequence ( MS) star 
(|Jura fe Helfandl [l98l IKastner fe Soken l200l . Deter- 
mining whether Mira B is a WD or a MS star is im- 
perative because its nature and evolutionary status im- 
pact our understanding of: 1) the efficiency of wind- fed 
accretion; 2) mass transfer through disks with sizes of 
~ 10 13 cm; 3) X-ray emission from symbiotic stars; and 
4) the generation of the bipolar structure in planetary 
nebulae (PNe). 

Because Mira AB is s patially resolved in t h e op- 
tical, X-ray, and radio dKarovska et all 119971 120051: 
iMatthews fe Karov ska 2006), it is an excellent place to 
investigate the accretion of a stellar wind from a de- 
tached companion. Roche lobe overflow is thermally un- 
stable when the donor star in a binary is an evolved 
giant with a mass in excess of 5M/6 (where M is the 
mass of the accretor; iWebbink. Rappaport. fe Savoniid 
1983). Hence symb iotic binaries, whic h typically have 
low- mass accretors (Mikolaicwska 2003), tend to trans- 
fer material via gravitational capture of the red-giant 
wind. In Mira AB, X-ray observations revealed a 
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bridge of material between the two stars (fkarovsk a et al.l 
2005) that likely arises from the gravitational effects, 
as the wind speed (6 km s _1 ; iKnapp fc Morrisl 119851 ) 
is of the same order as the stellar velocities. Mira 
may thus be the first clear example of a new mode 
of m ass transfer — dubbed 'wind R oche lobe over- 
flow' dMohamed fc Podsiadlowskill2007D or 'gravitatio nal 
focusing' (jde Val-Borro. Karovska. fc Sasselovl 120091 ) - 
that lies somewhere between standard Roche-lobe over- 
flow and Bondi-Hoyle accretion. This type of mass trans- 
fer may occur in those wide binaries tha t ultimately 
explo de as type la supernovae (SNIa; e.g.. iPatat et al.l 
l2007t iSimon et al.l [200l iBlondin et atll2009f ). with the 
yield from this channel of SNIa production depending 
on the efficiency of the mass transfer from the RG to 
the WD. Since the inferred accretion rate onto Mira B 
(M = LR/GM, where L is the luminosity due to ac- 
cretion, R is the radius of the accretor, and G is the 
gravitational constant) differs by almost two orders of 
magnitude depending upon whether Mira B is a WD or 
a MS star, determining the nature of Mira B is critical for 
estimating M and the ability of this channel to produce 
a SNIa. 

Mira AB will also test the hypothesis that the shapes 
of PNe are the result of interactions between a star 
in it s final stages of evolution and a binary compan- 
ion ()De Marcol I2009D . Observations indicate that pre- 
PNe commonly have bipolar morpholog ies and jets (e.g., 
IBuiarrabal et al.1120011: ISahai et aLll2007n . and most mod- 
els for prod ucing these structures require a binary com- 
panio n (e.g. JDe Marcd2009tlHuggins. Mauron. fc Wirthl 
120091 and references therein). iMiszalski et aLT ("2009). 
however, estimated that only 10-20% of PNe contain bi- 
nary central stars with separations small enough for the 
binary to hav e experienced a phase of common-e nvelope 
evolution, and lHuggins. Mauron. fc Wirthl ()2009f ) argued 
that the roundness of the majority of the halos of pre- 
PNe indicates that any binary companions in these sys- 
tems must either be quite distant (as in Mira) or have 
very low masses. Mira A is near its evolutionary end 
state, it has a distant binary companion, and recent 
ultraviolet (UV) images from the GALEX satellite re- 
yealed bipolar "stre ams" of knots to the north and south 
(|Martin et al.l 12007). Therefore, resolution of the issue 
of the nature of the companion would clear the path 
for models of Mira AB to confront questions such as 
how wide binaries generate bipolar structure in PNe, 
which stellar component launches the bipolar outflows, 
and when asymmetry first appears. 

In this paper, we use fast i?-band photometry to show 
that Mira B is a WD. We describe our rapid photomet- 
ric observations in Sj2j and the resulting light curves and 
power density spectra in S}3] In ijH we explain how the 
amplitude of the optical variations on time scales of min- 
utes reveal Mira B to be a WD. We also estimate a char- 
acteristic accretion rate onto Mira B and reinterpret the 
X-ray emission. We summarize our conclusions and ex- 
plore several implications in Sj5] 

2. OBSERVATIONS AND DATA REDUCTION 

A s part of a photometric survey o f 35 symbiotic stars 
(see lSokoloski. Bildsten. fc Holl2001l hereafter SBH), we 
observed Mira AB with the 1-m Nickel telescope at 
UCO/Lick Observatory on Mt. Hamilton, near San Jose, 



TABLE 1 
Log of Mira AB Observations 



Date, 


Obs. Start 


Duration 


t e 


At 


Num. 


U.T. 


(MJD) 


(hp) 


(s) 


(s) 


Pts 


1997 Sep 2 


50693.350 


3.6 


25.0 


46.996 


256 


1997 Nov 3 


50755.410 


1.6 


50.0 


70.0 


81 


1998 Aug 19 


51044.482 


1.0 


18.0 


38.990 


97 


1998 Aug 20 


51045.381 


3.5 


30.0 


49.988 


248 


1998 Sep 17 


51073.352 


4.6 


25.0 


46.0 


355 



Note. — Obs. Start is the start time for a given observation. 
Duration is the length of the observation. Columns 4 and 5 lists 
the exposure time and the time between exposure starts as t e 
and At, respectively. 



CA, on 5 nights between 1997 September and 1998 
September. Table [T] contains a log of the observations, 
each of which consisted of a series of exposures with an 
unthinned, 2048 x 2048 Loral CCD (referred to at Lick 
Observatory as 'CCD 2') that had 15-/xm pixels and a 
6'.3 x 6'.3 field of view. We chose exposure times, t e , of 
18 to 50 s to maximize the signal-to- noise ratio (S/N) 
while avoiding saturation of the program star or the one 
bright comparison star in the field of view (HD 14411). 
Keeping the exposure times at or above 18 s also kept the 
observing efficiency (t e /At, where At is the time between 
observation starts) from dropping much below ~ 50% 
and minimized scintillation noise, which can dominate 
the error at very short exposure times (SBH). Includ- 
ing chip pre-processing and readout times of between 20 
and 22 seconds, At ranged from 39 to 70 seconds. To 
obtain evenly spaced data points within a given obser- 
vation, and thereby produce data that were compatible 
with standard fast Fourier transform (FFT) routines, we 
employed a timing system developed especially for this 
project by W. Deitch (UCO/Lick). The resulting light 
curves spanned from 1.1 to 4.6 hours. To minimize the 
contribution of the red giant while still obtaining ade- 
quate flux from the bluer hot component to generate a 
high S/N on time scales of minutes, we used a Johnson B 
filter and observed only near the minima of the 332-day 
Mira pulse period, as displayed in Figure 1. 

For each observation, we reduced the CCD images us- 
ing standard techniques and performed aperture photom- 
etry, as described in SBH. For each image, our reduction 
included subtracting the electronics zero point and bias 
pattern and dividing by an average of multiple images 
of the sk y, using IDL so f tware based u pon IRAF rou- 
tines (see ICTillilandl [19921 : IMassevI [1991 . We produced 
differential light curves with respect to either the one 
bright comparison star in the field or a weighted, ensem- 
ble average of 2 - 4 comparison stars (depending upon 
whether the additional comparison stars improved the 
S/N enough to warrant the additional lost points due 
to cosmic rays). We estimated the background for each 
star from an annulus around that star. We omitted data 
points contaminated by radiation events ('cosmic rays'), 
points with extremely high background, and data taken 
during exceptionally poor weather conditions. None of 
the comparison stars showed any evidence for intra-night 
variability. 

3. ANALYSIS AND RESULTS 
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Fig. 1. — Long term visual light curve, from 1996 to 1999, from the American Association of Variable Star Observers (AAVSO), with the 
dates of our high-time-resolution B-band observations marked with arrows. Our observations all occurred near a photometric minimum of 
the pulsating red giant. Since we performed our third and fourth observations on subsequent nights (1998 August 19 and 20), they appear 
as a single arrow (the third from the left). 



TABLE 2 
Optical Variability of Mira AB 



Date 






Powerlaw 


Freq Range 




(U.T.) 


(mmag) 


(mmag) 


Index 


(mHz) 




1997 Sep 2 


2.7 


50.9 


-1.8 ±0.2 


0.08 - 3.6 


0.8 


1997 Nov 3 


2.6 


37.9 


-2.4 ±0.3 


0.2 - 3.4 


0.8 


1998 Aug 19 


2.6 


29.3 


-2.1 ±0.2 


0.3 - 9.3 


0.7 


1998 Aug 20 


2.3 


55.3 


-2.0 ±0.1 


0.08 - 5.5 


0.9 


1998 Sep 17 


2.5 


28.7 


-1.9±0.1 


0.06 - 10.8 


0.9 



a The measured rms variation, s, includes variations at all frequencies 
to which a given observation was sensitive; the range of frequencies 
sampled for each observation is approximately equal to the range 
over which we fitted the powerlaw model, shown in column 5. 

Rapid aperiodic variability, like the variability from 
cataclysmi c variables that is often termed 'flickering' 
(CVs; e.g.. lWarnerill995HBruclll2000r ). is evident on all 5 
nights. Figure [5] shows the differential light curves from 
the 5 observations. On each night, the measured root 
mean square (rms) variation, s, exceeded the rms varia- 
tion expected from the known errors, s oxp , by more than 
an order of magnitude. Table [5] lists s and s eX p (the es- 
timation of which is discussed in detail in SBH) for each 
observation. Figure [5] demonstrates that the £?-band flux 
can change by more than 5% (« 50 mmag) in less than 
5 minutes and more than 25% («0.25 mag) in less than 
2 hours. Taking into account the non- variable light from 
the red giant, which contributes « 50% of the flux at B- 
band (since Mira A and B had similar brightness in the 
HST F5 01N filter when Mira A was near pulse minimum 
in 1995. iKarovska et al.l 119971) , these rapid fluctuations 
correspond to variations of approximately 10 - 50% in 
the B-band light from Mira B. In summary, correcting 
for the contribution from the red giant, the i?-band flux 
from Mira B varied as fast as a few per cent (or tens of 
mmag) per minute, and by order unity in tens of min- 
utes, as exhibited in the light curve from 1998 August 20 
(Figure ©. 

The PDS of all 5 light curves reveal broadband power 
with a powerlaw shape — P v oc v a (where P v is the 



power density and v is the frequency) — between ap- 
proximately 0.1 and 10 mHz (see Table [2]). Since our 
data points were evenly spaced, we used the IDL rou- 
tine fft to compute the PDS. We normalized the PDS to 
the varianc e of the light curves usi ng the Miyamoto nor- 
malization (jMivamoto et al.lll991| ). in which the square 
root of the integrated power between two frequencies is 
equal to the fractional rms variation between those two 
frequencies. For PDS with a powerlaw index a = —2, 
the square root of the power at a given frequency equals 
the fractional rms variation summed over all frequencies 
above that frequency. Figure |4] shows the PDS corre- 
sponding to the light curves in Figure [2j Table [2] lists 
the powerlaw indices for each observation; they ranged 
from -1.8 to -2.4. The weighted average of the 5 pow- 
erlaw indices is a — —1.969 ± 0.004. Since this value is 
very close to a — —2, we refer to the stochastic variations 
from Mira as red noise. 

Extrapolating the PDS where necessary using a power- 
law index of a = —2, and taking into account the con- 
stant light from the red giant, we find that the average 
rms variation from Mira B between 0.1 and 10 mHz is 
approximately 9%. We derived this value from a mea- 
sured rms variation for the total -B-band light in this 
same frequency range of 4.5%, since emission from the 
red giant, which does not vary at these frequencies, con- 
stitutes approximately half of the B-band flux. Thus, the 
B-band rms variation from Mira B is at least twice the 
rms variation of the total B-band light. In the restricted 
frequency range of 1 - 10 mHz, which is useful for com- 
parison with published rms variations in this frequency 
range for CVs, the average rms variation from Mira B is 
3% (correcting a measured rms variation of 1.5% in this 
frequency range for the contribution from the red giant). 

We did not find any evidence for coherent, periodic 
variations on time scales of minutes to hours in the 
optical brightness of Mira B. Taking the broadband 
power into account (as d iscussed in the Appendix of 
iSokoloski fc Kenvon|[2003l ). no peaks in the PDS rose sig- 
nificantly above the level of the red noise. 
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Fig. 2. — Differential B-band light curves for Mira AB. The times between data points for the light curves are At =47, 70, 39, 50, and 
46 s, respectively. The size of the error bars is comparable to that of the point symbols. The time is U.T. on each date. 



4. DISCUSSION 

The nature of the accreting companion to Mira A may 
hold the key to improved understanding of mass transfer 
via a wind, as well as the way in which binaries gener- 
ate asymmetrical structure in PNe. Whereas e arly ob- 
servat ions of rapid CV-lik e optical variability by I Walker 
(|1957l ) and I Warner] (|1972f ) originally led most authors to 
assume that Mira B was a WD, later revelations that 
the X-ray luminosity w as lower than expected motivated 
Uura & Helfandl (fl98l and then lKastner fc Sokerl ([20041 ) 
to propo se that Mira B was in stead an accreting MS star. 
Though llreland et all (|2007l ) reported some marginally 
significant (4-er) features in the optical spectrum that 
appeared to suggest a MS compan ion, the UV spectrum 
is mu ch more indicative of a WD (jReimers fc Cassatellal 
1985). Below, we use a quantitative examination of the 
amplitude of the optical variability on time scales of 



minutes, along with an estimate of the expected opti- 
cal depth of the accretion-disk boundary layer, to show 
that the optical to X-ray observations are more consis- 
tent with a WD accretor. 

4.1. Rapid optical variations and the nature of Mira B 

The persistent presence of minute-time-scale stochas- 
tic optical variations with the observed amplitude is a 
strong indicator that Mira B is a WD. Qualitatively, 
whereas we are not aware of a single example of an 
accreting MS star with a high time resolution optical 
light curve resembling those of Mira, most accreting 
WDs have optical light curves that show Mira-like ape- 
riodic variations with a time scale of minutes — i.e., 
faster than the dynamical time of a MS star. Quanti- 
tatively, the PDS of our Mira light curves have shapes 
and strengths consistent with those of WDs and not 
MS stars. The PDS from Mira has a powerlaw index 
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Fig. 3. — Close-up view of the differential B-band light curve 
for Mira AB on 1998 August 20. The size of the error bars is 
comparable to that of the point symbols. The time is U.T. 

of approximately —2, like the high-frequency PDS of 
other accreting objects (including WDs, X-ray binaries, 
and AGN; e.g..lBruchlll992tlDone. Gierlinski, 'fc Kubotal 
[20071: iGierEnski, Nikolaiuk. fe Czernvl 120081) . We there- 
fore take the rapid variability to be related to accretion. 
Furthermore, the rms amplitude of the _B-band bright- 
ness variations from Mira B between 1 and 10 mHz (af- 
ter accounting for the contribution from the red giant) is 
around 3% — almost identic al to the characteristic value 
for CVs, which Barros (2008) found to have rms variabil- 
ities in the same frequency range of 2.4% in the g' and 
r 1 bands and 4.5% in the v! band. 

Theoretical considerations suggest that we should not 
see such rapid variations from a MS accretor. The dy- 
namical and viscous times at the inner edge of an ac- 
cretion disk around a MS star are » 100 times longer 
than those at the inner edge of a disk around a WD. 
Since these are likely the relevant time scales for setting 
the speed of brightness fluctuations, we expect an ac- 
creting MS star to exhibit stochastic variations that are 
stronger at low frequencies and weaker at high frequen- 
cies than the brightness fluctuations from an accreting 
WD. In fact, light variations have been detected from 
the accreting MS stars in Algol binaries, and these vari- 
ations prefere ntially occur on time scales of days rather 
than minutes; lOlson fc Etzell (|1995l ) found variations in 
the strength of double-peaked Ha emission, which they 
argue trace changes in the inner accretion disk, on time 
scales of days in a sample of 9 Algol-type binaries. Al- 
though we have no empirical constraints on variations 
from a MS-star disk in the frequency range of 1 - 10 mHz, 
we expect the amplitude of fluctuations in this frequency 
range to be orders of magnitude smaller than those from 
accreting WDs, if, for example, the PDS from an ac- 
creting MS star is simply that of a WD but shifted to 
lower frequencies. Therefore, given that the amplitude of 
the high-frequency fluctuations from Mira are consistent 
with those from CVs, and inconsistent with expectations 
for a MS star, we interpret the accreting object in Mira 
as a WD. 



4.2. Reinterpretation of the X-ray emission 

Before addressing the X-ray emission, we use the 
luminosity of Mira B in the UV and optical to estimate 
a characteristic accretion rate, M, onto Mira B. UV 
observations with the IUE satellite on 1980 June 16 
(when Mira A was at pulse minimum and therefore 
did not contribute to the UV flux) revealed a UV 
continuum luminosity between 1800 and 3200A of 
Luv ~ 4 x 10 32 (d/107pc) 2 erg s~ * (which we estimated 
from the top panel of Fig. 1 of IReimers fc Cassatellal 
119851 ) . As noted bv IReimers fc Cassatellal (|1985f ). this 
estimate is a lower limit to the UV luminosity, since an 
accretion disk could also radiate in the Lyman contin- 
uum, which IUE would not detect. Th e HST spectra 
of Mi ra A and Mira B in Figure 2 of fkarovsk a et al.l 
(|1997D . taken on 1995 December 11, shows that Mira B 
produced L opt « 4 x 10 32 (d/107pc) 2 erg s -1 between 
roughly 3200 and 5000A. Given that Mira B also emits 
at wavelengths longer than 5000A, we estimate that the 
total, accretion-powered luminosity of Mira B is typically 
at least L ps 10 33 (ri/107pc) 2 erg s -1 . This luminosity 
implies that Mira B accretes at a rate of at least J1%d ~ 
lO^Moyr- 1 (L/10 33 ergs" 1 ) (i?/10 9 cm) (M/O.6M )~ 

Consistent with our contention that Mira B is a WD, 
an accreting WD can naturally produce a UV spec- 
trum like the one seen from Mira. The IUE spec- 
tra from Mira between 1979 and 1983 contained lines 
with a range of ionization states, including N V 1240A, 
C IV 1550A, and Si IV 1400A, and line widths - 
from 135 km s" 1 for Fe II absorption, to few hundred 
km s" 1 for C IV 1550A to ~ 1 000 km s" 1 for Lya 
and Mg II (IReimers fc Cassatellal Il985l). Moreover, the 
Ly q profiles that iWood. Karovska. fc Raymond! (|2002l ) 
an d lWood fc Karovskal (|2004fl~ reconstructed from H2 flu- 
orescence lines observed by HST/ ST IS on 1999 August 
2 and by FUSE on 2001 November 27 had full widths 
at half maximum (FWHM) of approximately 1,100 to 
1,200 km s _1 . Velocities of greater than 1,000 km s _1 
are larger than the fastest Kepler velocities in the vicinity 
of an accreting MS star, but comparable to the Kepler 
velocities near an accreting WD. In fact, line widths of 
~ 1, 000 km s _1 are often seen in the UV spectra of CVs, 
which generally emanate from accretion-di s k win ds (e.g., 
iMauchd H99ll iFroning. Long, fc Baptistal [2003h . Fur- 
thermore. IReimers fc Cassatellal <| 19851 ) determined that 
the temperature of the material producing the high- 
ionization resonance lines in the IUE spectrum of Mira 
(which they take to be the coronal region of the disk) was 
on the order of 10 5 K — very reasonable for material near 
an accreting WD. 

Regarding the X-ray emission from Mira, an X-ray lu- 
minosity that is several orders of magnitude below the 
optical luminosity from the accretion disk, as we see 
from Mira B, is exactly what one would expect if the 
boundary lay er (BL) around the accretin g WD is op- 
tically thick (jPatterson fc Raymond! fl985f ) or when the 
spreading layer (SL) on t he accreting WD does not ex- 
tend above the thin disk (jPiro fc Bildstenl [20041 ). For a 
slowly rotating WD, roughly half o f the energy released 
via accretion is radiated by the BL (|Shakura fc Sunvaevl 
[1971 iL^nden-Bell fc Prjngid"ll97l or in the latitudinal 
spreading of the material after it arrives on the WD 
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Fig. 4. — Power spectra of the light curves shown in Figure [2] The average PDS in the lower right panel was produced by padding the 
light curves with enough zeros to generate frequency spacings that were the same for each of the nightly PDS. 



(jPiro fc Bildstenl 120041 ) . If the material is shocked as 
it arrives, the the high Keplerian velocities near the 
surface of the WD produce hard X-rays that escape if 
the plasma is optically thin. With an accretion rate of 
Mwd ~ 10 _10 M Q yr~\ we would expect a luminos- 
ity of at least a few times 10 32 erg s _1 from the BL 
around Mira B. Estimates of the X-ray luminosity, how- 
ever, all indica te that L T , is no more than a few times 
10 29 erg s- 1 (l.Tura fc Helfandl fl98i IKastner fc Soken 
[200llKarovska et al.ll2005Tr Although a detailed investi- 
gation of the X-ray emission is beyond the scope of this 
paper, the low L x most likely reveals that the boundary 
(or spreading) layer is optically thick and therefore emit- 
ting primari ly in the far- or extreme-UV instea d of in 
the X-rays. iWheatlev. Mauche. fc Matteil (120031 ) found 
that for the WD in SS Cygni, the BL became optically 
thick when the accretion rate rose above 10 -10 M Q yr _1 , 



in accord with our hypothesis that the accretion rate 
onto Mira B is high enough to quench the X-ray emis- 
sion. Since the accretion rate at which the BL tran- 
sitions from opt ically thin to optically th ick increases 
with WD mass (Popham fc Narayanl 119951 ). the optical 
thickness of the BL around Mira B suggests that the 
mass of Mira B is lower than that of the WDs in, e.g., 
the symbiotic binaries T CrB, RT Cru, and SS73 17, 
which contain WDs with masses greater than 1.0 Mq 
and BLs that are optically thin at accretion rates as high 
as a few times 10~ 9 M^ yr" 1 dLuna fc Sokoloskil 120071 : 
iKennea et"aIll2009HEze. Luna, fc Smithll2010l ). 

The existence of other WDs that accrete from the wind 
of a red giant and produce a very low X-ray luminos- 
ity supports our claim that the X-ray luminosity from 
Mira B is consistent with WD accretion. For exam- 
ple, the accretor in R Aqr, which Jura fc Helfandl 0984) 
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also suggested was a MS star based on its low X-ray 
luminosity, was later shown by IU E to have an effec- 
tive temperature of 61 ,000 K and the ra dius of a WD 
( Meier fc Kafatod fl995T ) . More generally. iMiirset et alJ 
(1997) found that the bolometric luminosities of the 
wind-fed WDs in symbiotic stars (such as EG An- 
dromedae, PU Vulpecula, and RX Puppis) often exceed 
their X-ray luminosities by as much as 5 orders of mag- 
nitude. 

Finally, a different suggestion for Mira B — that 
the faint X-rays are coronal emission from a mag- 
netically active, rapid ly rotating low-mass MS star 
(|Kastner fe Sok cr 2004) — cannot explain either the 
strong UV emission (compared to X-rays) or the red- 
noise-type optical variations. Whereas the luminosity of 
the rapidly variable component of the optical emission 
from Mira B is on the order of 10 32 erg s _1 , the time- 
averaged power of optical r adiation from stella r flares 
rarely exceeds 10 30 erg s _1 (Shakhovskaval [T989l ). More 
importantly, a MS star producing an X-ray luminosity of 
a few times 10 29 erg s _1 due to coronal emission would 
only have a tim e-averaged optical p ower from flares of 
~ 10 29 erg s" 1 (jShakhovskaval H989) ■ We can thus rule 
out coronal emission from a magnetically active, rapidly 
rotating, MS star on energetic grounds. The radio emis- 
sion from Mira B is also one to three orders of magnitude 
greater than expected from a magnetically active dwarf 
([Matthews fc Karovskal 120061 ) . In contrast, all available 
observations are consistent with the picture in which a 
WD is accreting from the wind of the red giant. 

5. CONCLUSIONS AND IMPLICATIONS 

Our main conclusion is that the rapid optical variabil- 
ity from Mira AB reveals Mira B to be a WD. The ampli- 
tude of the optical brightness fluctuations from Mira B on 
time scales of minutes is the same as those from accret- 
ing WDs in CVs, and significantly larger than one would 
expect from an accreting MS star. The UV spectra — 
especially the ionization states and widths of the lines — 
are consistent with the conclusion that Mira B is a WD. 
The X-ray luminosity and spectrum are also consistent 
with this conclusion if the accretion-disk boundary layer 
is optically thick, as one would reasonably expect for our 
estimated accretion rate of ~ 10~ 10 M w yr -1 (e.g., see 
iWheatlev. Mauche. fc Matteill2003[ ). 

This accretion rate, which we infer from the opti- 
cal and UV luminosities of Mira B, is consistent with 
constraints on the effective temperature of the WD 
from UV spectra. Sion (1999) noted that prolonged 
accretion will re-heat an otherwise cooling WD, mak- 
ing it potentially visible in the UV or optical contin- 
uum. Such rejuvenated WDs are seen in dwarf novae 
systems during quiescence, reveal ing effective tempera- 
ture of T eff = 10.000-20.00 K dTownslev fc Gansickel 
20091). This WD re-heating (|Townslev fc Bildstenl 12003 
Townslev fc Gansickel [2009) and the resulting T e s de- 
pend on both M and the time-averaged accretion rate, 
(M), thereby providing another constraint for our anal- 
ysis of Mira B. From the UV continuum flux upper limit 
Fa(1250A) < 10~ 14 ergs cm -2 s _1 , Reimers and Cas- 
satella (1985) found T eS < 14, 000 K for R = 0.012i? Q 
(for a WD with M = 0.6M Q and a distance of 77 pc; 
a slightly higher T e g is allowed at our adopted distance 
of 107 pc). A direct application of iTownslev fc Bildstenl 



(2003, see their Figure 1) converts this upper limit on 
T e ff into the long-term accretion rate constraint (M) < 
2 x 10- 10 Af Q yr _1 , in agreement with our earlier esti- 
mate. 

At this accretion rate (M « 10- 10 M Q yr -1 ), the ac- 
cumulated matter on the WD will burn unstably and 
cause a classical novae, leading to the violent ejection 
of ~ 1O" 4 M of material at > 1000 km s _1 . This 
ejected mass is larger than the total mass of the wind 
between the WD and the red giant, leading to little de- 
celeration of the shock as it plows through the wind in 
the few months following the eruption. A nova erup- 
tion from Mira would thus perhaps look more like those 
of the symbiotic slow novae than, say, an X-ray bright 
eruption of the symb iotic recurrent nova RS Ophiuchi 
(|Sokoloski et al.ll2006[ ) or the 7-ray producing ex plosion 
of the symb iotic V407 Cvgni (lAbdo et alj |20ldh. From 
the work of ITownslev fc Bildstenl (|2004l ). we would ex- 
pect a nova recurrence time of « 10(2) x 10 5 yr for 
M = 0.6(1. 2)M Q and (M) = 2 x 10- 10 Af o yr" 1 . 

Our conclusion that Mira B is a WD accreting 
at M ~ 10- 10 A/q yr -1 has implications for the 
mode and effic iency of mass transfer in wide bina- 
ries. Although |Mohamed fc Podsiadlow ski (2007}) and 
Ide Val-Borro. Karovska. fc Sasselovl (|2009l ) have sug- 
gested that a focused wind could enable up to tens 
of percent of the donor wind to be captured, the 
accretion rate needed to produce Mira B's luminos- 
ity and effective temp erature is just 0. 1% of the RG 
wind mass- loss rate (jRyde et alJ [2000). It is also 
below the Bondi-Hoyle accretion rate of Mbh ~ 
nr 8 M yr- 1 (A//0.6M Q ) 2 (7kms- 1 /-y oo ) 3 , where is 
the relative velocity of the wind at Mira B. Thus we find 
no evidence that either gravitational focusing or wind 
Roche-lobe overflow is enhancing the accretion efficiency 
onto the WD in this binary, despite the image showing a 
bridge of X-ray emitti ng material between the two stars 
(|Karovska et all 120051 ). 

Identifying Mira B as a WD also has ramifications 
for the launch site of the bipolar "streams" that have 
been imaged in t he near- and far-UV with Galex 
( Martin et al J 120071) a nd in the optical from the ground 
( Meaburn et al.1 20091 ) . Based on the behavior of other 
WDs that accrete from wind-fed disks in symbiotic bina- 
ries, we speculate that the the bipolar streams emanate 
from the accreting W D rather than, or in a ddition to, the 
pulsating AGB star. iMikolaicwska (2003) has suggested 
that symbiotic stars in which the wind of the red giant 
is focused toward or bital plane — as has been seen to 
be the case for Mira (|Karovska et al.|[2005D — preferen- 
tially contain large, unstable accretion disks, which ten d 
to produce sporadic bipolar outflows (jSokoloskil [2003) . 
Galex observations indicate that the streams from Mira 
contain about 10~ 7 M p of material (Don N eill, private 
communication), and iMeaburn et al.1 (|2009l ) estimated 
that the streams were ejected over the course of the past 
~ 1,000 yr. The average rate of flow into the streams 
was thus roughly (M str cams) ~ 10- 10 Af Q yr -1 . Given 
our estimate of the accretion rate, Mira B could fairly 
easily pump ~ 10~ n Mq yr -1 into the streams, which is 
roughly consistent with t he observational estimate of the 
mass-loss from the WD bv lWood. Karovska. fc Raymond! 
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( 2002) during the WE era. Given the large uncertainties 
on both the rate of accretion onto Mira B, the average 
rate of mass loss from Mira B, and the rate of flow in 
the bipolar streams, we conclude that the WD remains a 
viable the launch site, and that the streams do not nec- 
essarily originate from the RG. Mira AB therefore does 
not necessarily require a mechanism for generating fast, 
bipolar outflows from the AGB star itself, and it does 
not provide evidence that such mechanisms are needed 
to generate bipolar symmetry in PNe. 
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